Magnetic anisotropies of late transition metal atomic clusters 
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We analyze the impact of the magnetic anisotropy on the geometric structure and magnetic or- 
dering of small atomic clusters of palladium, iridium, platinum and gold, using Density Functional 
Theory. Our results highlight the absolute need to include self-consistently the spin orbit interac- 
tion in any simulation of the magnetic properties of small atomic clusters, and a complete lack of 
universality in the magnetic anisotropy of small-sized atomic clusters. 
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Nanostructures of all kinds display a wealth of fascinat- 
ing geometric, mechanical, electronic, magnetic or optical 
properties. The exploding field of Nanoscience pretends 
to understand, handle and tailor these properties for hu- 
man benefit. Atomic clusters and chains, molecular mag- 
nets, and a number of organic molecules like for instance 
mctallocencs indeed show novel magnetic behaviors, that 
include the enhancement of magnetic moments due to the 
reduced coordination and symmetry of the geometry [1J, 
and a rich variety of new non-collinear magnetic struc- 
tures that are absent in bulk materials [2J. Among all 
these devices, metallic atomic clusters (MACs) [3] stand 
out since, on the one hand they represent the natural 
bridge between atomic and materials physics and, in the 
other, they can be grown, deposited on surfaces or em- 
bedded in diverse matrices, and characterized with rela- 
tively well-established techniques. Further, the magnetic 
properties of MACs show promise for a wide spectrum of 
applications, ranging from medicine to spintronics. 

Magnetism in small MACs has been extensively stud- 
ied both experimentally and theoretically along the past 
decade |U |5]. It is therefore somewhat surprising that, 
even though spin-orbit effects are also expected to be en- 
hanced in MACs, there are very few published theoreti- 
cal papers that include the spin-orbit interaction (SOI) in 
their simulations [BJ |5J H] ■ Among these, only a handful 
have actually looked at the explicit effects of the SOI on 
the magnetism of MAC. Pastor and coworkers have stud- 
ied the magnetic anisotropy (MA) [10] of clusters of 3d 
Transition Metal atoms, using a phenomenological tight- 
binding scheme. To the best of our knowledge, there 
are no ab initio studies of the impact of the SOI on the 
magnetism of 4d and 5d MAC. This includes not only 
explicit calculations of the magnetic anisotropy energy 
(MAE) , but also and more importantly, whether the SOI 
modifies the ground state and the tower of lowest lying 
(magnetic) states of a given cluster. 

We present in this article a series of calculations of se- 
lected atomic clusters of 4d- and specially 5d-elements, 
that show that the SOI is a key ingredient in any ab 
initio simulation of heavy transition metal MAC, that 
should not be overlooked. These include Pd„, Ir„, Pt„ 
and Au„, with n=2, 3, 4 and 5, and also Pdg and Aug, 



AU7. The rationale behind our choice is that 5d MAC 
are expected to have the largest MAEs; among all 5d 
elements, gold is monovalent and should display a sim- 
pler behavior; comparing the magnetic properties of gold 
MAC with those of its neighboring Pt and Ir elements 
should shed light on the possible trends; further, Pd is 
isoclcctronic to Pt and it is important to check whether 
they display similar magnetic behavior. We have found 
a number of interesting results. First, the SOI always 
alters slightly the geometry of MAC, therefore generaliz- 
ing the results obtained for platinum[9j to other elements. 
Furthermore, the SOI can lead to qualitative changes in 
the atomic geometry in specific cases, as for example Ir3, 
Ir 4 and AU7. Second, the MA generated by the SOI has 
a strong impact on the magnetism of the tower of lowest 
lying states of the clusters. The MA selects orientations 
that are typically parallel (PR) or perpendicular (PP) 
to specific bonds in each cluster. Thus, MA and ten- 
dency to ferromagnetic alignment conspire together for 
ferromagnetism, or frustrate each other, depending on 
the geometrical arrangements of the atoms in each clus- 
ter. Hence, the lack of universality of small MACS is 
even more evident for their magnetic properties than for 
their geometry. Finally, atomic spin and orbital moments 
align in parallel to each other. 

We have used the ab initio Troullier-MartinsfTT] pseu- 
dopotential code SIESTA [12]. We have added to this 
code a fully self-consistent and spin non-collinear imple- 
mentation of the SOI as described in detail in Ref. |13j . 
This implementation has been tested satisfactorily in sev- 
eral bulk semiconductors and metals [13 a , some molecular 
magnets [14] . a large number of atomic chains[15 and a 
good number of small metallic (in the present work) and 
intermetallic clusters [13] by now. We have crosschecked 
our results for Ir2 and Pt2,3 with the plane-wave, Van- 
derbilt pseudopotentials code Quantum ESPRESSO [16]. 
We have employed both LDA and GGA approximations 
[T71 1TB] for all dimers and trimers to assess the reliabil- 
ity of our results, except for Ir 3 . We have used LDA 
for larger MACs since LDA provides a slightly more ac- 
curate description of bulk Pd, Ir, Pt and Au, includ- 
ing the magnetism of Pd[19]. We proved recently that 
the ground state and first excited isomers of Pd clus- 



Typeset by REVTeX 



2 



ters are rather insensitive to the choice of correlation 
functional and pseudopotential[20]. We have found here 
that the MAE further enhances the consistency of LDA 
and GGA in Pd, Pt and Ir MACs. In contrast, we con- 
firm that the choice of functional is important for Au„ 
clusters [5J [H], even after including the SOI. We have bor- 
rowed the atomic configuration for the pseudopotentials 
from Quantum ESPRESSO package's web page [16]. We 
have checked that they produce a correct description for 
the bulk phases, and have also benchmarked with other 
choices of matching radii in the case of some Pd„, Pt„ 
and Ir„ MACs. We have included non-linear core correc- 
tions for Pd to account for the effect of semicore states 
in the valence. We have used in all cases a very com- 
plete triple zeta doubly polarized basis set (three pseudo 
atomic orbital for s and d shells, and two for the unoc- 
cupied p shell), with long radii that provides a very ac- 
curate description of the structural and electronic prop- 
erties of the corresponding bulk materials. We have ex- 
plicitly tested and set tolerances in all parameters (den- 
sity matrix, electronic temperature, real space grid and 
forces) to achieve accuracies in the range 1-3 K. We have 
calculated the individual atomic spin moments Si using 
a Mulliken analysis, while to estimate the atomic orbital 
moments L t we needed to compute previously the ma- 
trix elements of the orbital angular operator in the cho- 
sen pseudo atomic basis. We have taken between two 
and five different seeds both for the geometry and for the 
magnetism of each cluster, making a total of five to ten 
different initial configurations for each of the simulated 
clusters. 

We discuss first the highly non-trivial behavior of the 
magnetic anisotropy of the dimers Pd2, Ir2 and Pt2 and 
Au2, that we show in Figs. 1 (al) through (dl) for LDA. 
We have found that GGA agrees with LDA in all the 
essential details 22J. Pd2 shows a weak MAE of about 
5 meV that favors an alignment perpendicular to the 
dimer axis at its equilibrium distance. This dimer dis- 
plays a low to high orbital angular momentum transition 
at longer elongations, which swaps the MA to PR. In 
marked contrast to Pd2, the spin moments of the isoelec- 
tronic Pt dimer show PR alignment at its equilibrium 
distance. The MAE is estimated as 100 meV or 35 meV, 
depending on whether LDA or GGA is used. The two 
electrons of a single Pd or Pt atom arrange in a high spin 
(S=l) configuration, according to Hund's rule. We find 
that the impact of the SOI is larger for the Pt2 dimer, 
where the spin moment is noticeably smaller than 2 (1b, 
and different for the two orientations, as shown in Fig. 1 
(c2). For L2, the SOI also favors PR alignment, but the 
predicted MAE at the equilibrium distance of the dimer 
is of about 60 meV in the LDA approximation. The 
spin moment is twice larger than for Pt2 and Pd2 and is 
also slightly different for the two orientations. We have 
found that the two electrons in the gold dimer pair in a 
singlet state and therefore this clusters does show zero 
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FIG. 1: (Color online) Calculated ground state energy E, spin 
moment per atom times the gyromagnetic ratio g Si and or- 
bital moment per atom Li of (al-a3) Pd2, (bl-t>3) Ir2, (cl-c3) 
Pt2 and (dl-d3) Au2, as a function of inter-atomic distance 
for spins perpendicular (solid line) or parallel (dotted line) 
to the axis of the molecule. A different reference energy has 
been set for each dimer to provide for uniformity in the en- 
ergy values. Figures (d2-d3) demonstrate that the moments 
of the gold dimer are exactly zero. 



MAE. Figs. 1 (a3-d3) show the orbital moments of the 
different dimers. The atomic orbital moments align par- 
allel to the spin moments always, therefore spelling out 
that the L ■ S product favors ferromagnetic alignment of 
these moments. Notice that the SOI enhances the orbital 
magnetism in these dimers, as is reflected by the fact that 
the orbital moment of Ir 2 and Pt 2 is similar, but much 
larger than that of Pd2- If we take these results as a 
rough guide for the anisotropy of bigger clusters, we ex- 
pect a tendency to PP alignment for Pd MACs, while Pt 
and Ir clusters will prefer PR orientations. Before mov- 
ing to bigger clusters, we note that a tendency to PP 
alignment is easily fulfilled for planar structures, while it 
leads to frustration (e.g. non-collinearity) if the geom- 
etry is three-dimensional; tendency to PR alignment on 
the other hand can only be fulfilled without frustration 
for linear geometries. 

We pass now to discuss the geometry and magnetism 
of the ground state of the trimer MACs, that we show 
in Fig. 2. We find that Pd and Pt MACs form a tri- 
angle, while Ir arranges in a straight line, and AU3 ar- 
ranges as either a triangle (LDA) or a zigzag (GGA). 
The tendency to PP alignment of the spin moments in 
Pd2 leads to expect out-of-planc alignment for the Pd 
triangle. But the moments in Pd3 actually align ferro- 
magnetically in the plane of the triangle, in a fork- like 
fashion. A closer look to this cluster reveals that the in- 
teratomic distances in the trimer are larger, and actually 
occur at the crossover length from PP to PR alignments 
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FIG. 2: Equilibrium geometry and spin moments g Si in the 
ground state of the trimers (a) Pd3, (b) h-3, (c) Pt3 and (d) 
Au 3 . 



in the dimer. This leads to a sort of magnetostrictivc 
effect whereby one of the bonds becomes slightly longer 
than the other two, hence slightly distorting the equilat- 
eral triangle. The spin moments align perpendicular to 
that bond. Pt3 is a perfect equilateral triangle, that does 
shows no magnetostriction at all. The spin moments in 
platinum arrange in a fork-like fashion in the plane of 
the triangle, as those in Pd3, but now the moments are 
slightly non-collinear. Pd3 and Pt3 also share the first 
excited isomer, which shows out-of-plane anisotropy and 
ferromagnetic alignment, while the second one is a vortex 
for Pd 3 , and shows out-of-plane antifcrromagnetic align- 
ment in the case of Pt3. A rough estimate of the MAE 
for these two trimers can be made by taking the energy 
difference between the ground state and the first and sec- 
ond excited isomers, since they share the same geometry. 
This criterion yields MAEs of 25 and 35 meV for Pd and 
Pt respectively, if we use the first state, and 60 and 90 
meV when we use the second state. Interestingly, these 
estimates show similar values for the MAE of these two 
MACs, despite their different atomic weight. The PR 
MA in Ir 3 is strong enough to modify its geometry, sta- 
bilizing a linear structure instead of the triangle that is 
found when the SOI is not included. The first excited 
state in L3 is a paramagnetic linear chain that has an 
excitation energy of 7 meV, while the second, that lies at 
a much higher energy of about 300 meV is a triangle. AU3 
is weakly magnetic. LDA predicts a triangle with clear 
PP in-plane anisotropy which leads to a non-collinear 
state, while GGA yields a zigzag arrangement, with the 
spin moments also showing PP in-plane anisotropy. The 
first excited state of LDA AU3 is a triangle that shows 
out-of-plane PP anisotropy and has a MAE of 8 meV. 
GGA predicts smaller MAEs for gold than LDA, of the 
order of a few Kelvin |22J. 

We confirm that the atoms in the ground and first 
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FIG. 3: Equilibrium geometry and spin moments g Si in the 
ground state of the tetramers (a) Pd4, (b) Ir4, (c) Pt4 and 

(d) Au 4 . 



excited states of Pd4 arrange as a tetrahedron [20]. [22] . 
On the contrary, all lowest-lying states of the fifth-row 
tetramers show planar geometries, as shown in Fig. 3 
for the ground state clusters. Platinum tetramers form 
rhombuses, while the lowest lying Ir 4 isomers are per- 
fect squares. An exception is the second excited state 
of iridium, that is slightly bent in a sort of saddle. Gold 
displays a richer variety of lowest lying planar structures. 
The tendency to both ferromagnetism and PP MAE in 
the Pd 4 MACs leads to strongly frustrated non-collinear 
states, with small energy differences among them. The 
ground state shows the magnetostrictive effect found in 
Pd3, while the magnetic arrangements in the excited 
states precludes it. The ground state of the iridium 
and platinum tetramers shows in-plane MAE, where the 
atomic spin moments are aligned along one of the diag- 
onals of the square or rhombus. Their first excited state 
display out-of-plane MAEs, on the contrary. The excita- 
tion energies of these lowest lying states are of about 130 
and 50 meV for Ir 4 and Pt 4 , respectively. Ir 4 and Pt 4 
MACs do not show magnetostrictive effects. The lowest- 
lying states of gold clusters are paramagnetic, except for 
the second one, that shows a ferromagnetic alignment 
with out-of-plane MA. The excitation energy of the first 
gold excited state is of 260 meV[2"2"]. 

The tower of lowest lying states of the palladium, plat- 
inum and iridium pentamers show fully 3d arrangements. 
Palladium MACs typically form triangular bipyramids. 
The ground and first excited state of Ir5 are Keops pyra- 
mids, while the second excited state is planar. The 
ground state of Pts is also a Keops pyramid, but the two 
first excited states are triangular bipyramids [22]. Gold 
pentamers are all planar; the ground state and first ex- 
cited isomer displaying an arrangement reminiscent of 
the Olympics logo. The atomic spins in all Pds states 
are ferromagnetically aligned and display the same clear 
tendency to PP MAE as in Pd3 and Pd 4 . This is to be 
contrasted to the case of platinum and iridium MACs, 
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FIG. 4: Equilibrium geometry and spin moments g Si in the 
ground state of the pentamers (a) Pds, (b) Irs, (c) Ptg and 
(d) Au 5 . 

where the tendency to PR MAE can not be satisfied well 
in these 3d structures. This effect is specially acute for 
the Keops geometry. The trade-off in this case between 
spin ferromagnetism and PR alignments gives rise here 
to small non-collinearities, as shown in Fig. 4. The spin 
moments in these two MACs lie in the basement's plane. 
They point along its diagonal for Pts, while for Irs the 
moments are aligned parallel to the atomic bonds. The 
ground state of Aus is magnetic and displays in-plane 
MAE. The favored anisotropy is actually PP, which can 
not be satisfied, giving rise to a non-collinear arrange- 
ment, that we show in Fig. 4. The first excited state 
of Aus shows out-of-plane anisotropy, and has an excita- 
tion energy of only 4 Kelvin, revealing again the extreme 
softness of small gold MACs [22]. 

We comment finally on the behavior of a few MACs of 
larger size. The ground and first excited states of Pdg 
are square bipyramids. The SOI gives rise to small spin 
and angular moments in the ground state, that ceases to 
be paramagnetic. It is the first state, whose excitation 
energy is of 6 meV, which becomes paramagnetic. The 
ground state of Aug is degenerate. One of the states is a 
paramagnetic planar triangle while the other is a square 
bipyramid, that has small (non-collinear) spin moments, 
see Figs. 5 (a) and (b). Finally, we have found that the 
SOI stabilizes a planar atomic arrangement in Au 7 . The 
ground state presents two degenerate magnetic configu- 
rations, with in-plane and out-of-plane MAEs, that we 
show in Figs. 5 (c) and (d). These again spell out the 
extreme softness of gold clusters [22] . 

To summarize, we have performed a theoretical analy- 
sis of the impact of the spin-orbit interaction in the mag- 
netic properties of small clusters of Pd, Ir, Pt and Au. 
We have found that the SOI modifies, slightly always and 
qualitatively in a few cases, the geometry of these clus- 
ters. More importantly, the spin-orbit interaction gives 




FIG. 5: Equilibrium geometry and atomic spin moments g 
Si of the two degenerate ground states of Au6 ((a) and (b)) 
and Au7 ((c) and (d)). 



rise to magnetic anisotropics. These sometimes enhances 
the tendency to ferromagnetic ordering, while other com- 
pete with it, yielding some non-collinear structures. The 
MAE also gives rise to an interesting magnetostrictive 
effect in Palladium clusters. Our results highlight a lack 
of universality in the magnetic anisotropy of small-sized 
atomic clusters. 
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